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Abstract. We study the cosmological evolution based upon a £>-dimensional action 
in low-energy effective string theory in the presence of second-order curvature 
corrections and a modulus scalar field (dilaton or compactification modulus). A 
barotropic perfect fluid coupled to the scalar field is also allowed. Phase space analysis 
and the stability of asymptotic solutions are performed for a number of models which 
include (i) fixed scalar field, (ii) linear dilaton in string frame, and (Hi) logarithmic 
modulus in Einstein frame. We confront analytical solutions with observational 
constraints for deceleration parameter and show that Gauss-Bonnet gravity (with no 
matter fields) may not explain the current acceleration of the universe. We also study 
the future evolution of the universe using the GB parametrization and find that big 
rip singularities can be avoided even in the presence of a phantom fluid because of the 
balance between the fluid and curvature corrections. A non-minimal coupling between 
the fluid and the modulus field also opens up the interesting possibility to avoid big 
rip regardless of the details of the fluid equation of state. 
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1. Introduction 

During the past few years precision cosmology experiments have confirmed the big 
bang/inflationary standard scenario of a flat universe characterized by nearly scale- 
invariant, Gaussian primordial density perturbations. Two thirds of the total energy 
density, which is close to the critical one, is contributed by an exotic form of matter 
popularly known as dark energy (see [2| and references therein for an overview 
on the subject). Moreover, recent observations of extra-galactic supernovae treated 
as standardized candles indicate that the universe is presently undergoing a phase of 
accelerated expansion which receives an independent support from CMB and large-scale 
structure data. Within the framework of General Relativity, the cosmic speed-up can 
be fueled by an energy-momentum tensor with negative pressure - dark energy. The 
analysis of high redshift type la supernovae exhibits a marginal evidence of a possible 
super-acceleration j^J. The latter eventuality has triggered many speculations on the 
underlying nature of dark energy, among which the phantom hypothesis has been under 
active investigations (see Refs. [7| JO] for early works on phantom). By definition, 
a phantom is a fluid with pressure and energy density satisfying an equation of state 
p < —p; if the fluid is represented by a scalar field, its kinetic energy has a negative 
sign (hence phantoms and ghosts are relatives). This kind of contribution leads to 
a finite-time future singularity dubbed big rip. We also note that there exist other 
types of sudden future singularities even in the absence of the phantom [H| [IS]- ^n 
updated list of references on these subjects, including first and latest proposals, big rip 
phenomenology, and criticism, can be found in [HJ . 

The search for a viable semiclassical or quantum model for the dark energy 
component and its dynamical properties is one of the topics explored by "string 
cosmology", which attempts to understand the cosmological implications of string 
effective actions. In particular, quantum corrections to the Einstein-Hilbert action 
naturally arise for the closed string, bosonic or supersymmetric, and modify the 
gravitational interaction in a way which might be testable in the near future. The 
Gauss-Bonnet (GB) combination of curvature invariants with fixed dilaton is of special 
relevance: (z) it represents the unique 5D leading-order extension of the Einstein- 
Hilbert action that leads to second order gravitational field equations, and (ii) the 
highest derivative occurs in equations linearly, thereby ensuring a unique solution; finally 
(Hi) with the GB parametrization, graviton interactions are ghost-free and spacetime 
perturbations are wave-like. However, in four dimensions, the GB term with fixed 
dilaton is topological and does not contribute to the equations of motions. Then, the 
only way to modify the cosmological evolution is to consider GB term in the bulk in a 
braneworld setup, so that the induced 4D Friedmann equation on the brane depends on 
the bulk back-reaction. The GB braneworld and its observational consequences, both 
in the early and recent universe, are still under extensive study (see for a list of 
references) . 

Second-order gravity terms can generate non-trivial effects in a four- dimensional 
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spacetime even outside the braneworld lore, provided that we relax some common 
assumptions: the most important ones are fixed dilaton and GB parametrization. In 
the first case, it is assumed that non-perturbative potentials arise and dilaton and other 
moduli freeze out with a minimum non-zero mass, so that to preserve the observational 
bounds for the gravitational interaction PQEj- However, a massless dilaton still can 
be compatible with observations [T7J [TH] EH] (see also [2T)1 l2"T| ). 

If the dilaton varies in spacetime, gradient terms appear and the higher-order 
curvature invariants interact with a non-constant a'-coupling. Also, compactification 
of the 26 or 10-dimensional target space is encoded in residual modulus fields which in 
general will evolve in time. Both cases result in a dynamical field non-minimally coupled 
to gravity. 

Although the Gauss-Bonnet combination of curvature invariants is ghost-free, we 
should bear in mind that in the high-energy limit the uniqueness of the GB term is not 
guaranteed and, from the point of view of string theory, the GB parametrization cannot 
be distinguished from the others (see references in |22j). 

Dilaton and modulus cosmologies in higher-order gravity were considered, e.g., in 
|22j-|39j. As regards the recent universe, dark energy models with higher-order terms 
were inspected, e.g., in j^U] i n the case of fixed moduli, while a dynamical dilaton as a 
dark energy candidate was studied in jUlH^l- Recently Nojiri et al. [13] investigated a 
dark energy scenario in GB gravity coupled to a dynamical scalar field with non-negative 
potential. 

The understanding of the impact of higher-order curvature terms on the 
macrocosmic evolution may open up the possibility to verify high-energy fundamental 
theories through cosmology. We address this issue in two ways. First, we find 
cosmological solutions in the absence of an extra fluid in the action and constrain them 
via the recent supernovae data jH]. It turns out that, while solutions with non-GB 
parametrizations can fit the experimental evidence, the GB solutions are in general 
not viable, unless one takes the perfect fluid into account or, in the case of finite- 
time future singularity solutions, one fine-tunes the parameters of the model. Second, 
we solve numerically the equations of motion with a perfect fluid in four dimensions 
and study the future behaviour of its energy density and the Hubble parameter. In 
particular we are interested in the case of a phantom fluid which leads to the growth 
of the energy density of the universe. Then higher-order curvature corrections become 
inevitably important around the big rip. In fact, we will show that both second-order 
curvature terms and a non-minimal coupling between the fluid and the modulus field 
deeply affect the cosmological evolution. We shall classify the situations in which big 
rip singularities are reached or avoided. We keep arbitrary coefficients in the second- 
order gravity action, considering then the usual Gauss-Bonnet parametrization as a 
special case. Our results are in agreement with past investigations in literature where 
comparison is possible. Also, this work extends previous analyses to a spacetime of 
arbitrary dimension (including the 4D case) with a (coupled) barotropic fluid being 
taken into account. 
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The paper is organized as follows. In section |2] we present a general action which 
involves curvature, scalar field, and barotropic fluid terms, and derive basic equations 
for this system. 

In section 0] we take up the issue of de Sitter and inflationary solutions in the case 
of a constant scalar field. We then work out the phase space and stability analysis for 
a linearly changing dilaton in the string frame in section 0] Some technical material is 
confined to the Appendix. 

Section is devoted to the logarithmic modulus case with stabilized dilaton, which 
solves either exactly or asymptotically the equations of motion. The stability of these 
solutions is studied, together with their application as concrete models for the late time 
evolution of the universe. 

The properties of second-order gravity in the context of (phantom) dark energy are 
discussed in section El We show that it is possible to avoid big rip singularity in the 
presence of the Gauss-Bonnet corrections with a dynamically changing modulus field. 
The coupling Q between the modulus and the fluid can also play important roles in 
determining the future fate of the universe. Conclusions are given in section [7| 

2. General setup 

In our notation, is the D-dimensional metric with signature (— , +,•••,+) and Greek 
indices run from to d = D — 1. The Riemann and Ricci tensors are defined as 



R a »» = d„r% - d v V% u + r« Y% - r^r^ and R„ v = R« respectively, where Y% is 



the Christoffel symbol. The gravitational coupling kd and the Regge slope a' are set to 
unity. 

2.1. General action 

The gravity+matter action we begin from is (e.g., [4*4*] 143]) 



where g is the determinant of the D- dimensional metric, is a scalar field corresponding 
either to the dilaton or to another modulus, and / is a generic function of the scalar 
field and the Ricci scalar R. u, £ and V are functions of <fi. In this paper we do not 
consider the cosmological dynamics in the presence of the field potential V, but we 
include this term in deriving basic equations. C p is the Lagrangian of a .D-dimensional 
perfect fluid with energy density p and pressure p. Later on we shall assume that the 
barotropic index w = p/p is a constant. In general the fluid will be coupled to the scalar 
field 0. 

Finally, a'-order quantum corrections are encoded in the term 




d D x^g- R) - ^(0)(V0) 2 - V(0) + + 4 0) 



(1) 



Cf = a x R a ^ v R a ^ + a 2 R pv R^ + a 3 R 2 + a 4 (V0) 4 , 



(2) 
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where dj are coefficients depending on the string model one is considering. In section I2~31 
we will specify the functions /, u, V and £ together with the coefficients a« according 
to the scenarios of interest. 

We assume that the target spacetime is described by a flat Friedmann-Robertson- 
Walker metric 

d 



ds 2 = -N 2 {t)dt 2 + a 2 {t) ^(d^) 2 , 



(3) 



i=l 



where N(t) is the lapse function (N — 1 in synchronous gauge) and a(t) is the scale 
factor describing the physical size of the universe. The size of the causal connected region 
centered on the observer is given by the inverse of the Hubble parameter H = d/(Na). 
In the following, dots and primes denote derivatives with respect to synchronous time t 
and (ft, respectively. 



2.2. Equations of motion 

With the metric (JHJ), the Riemann invariants read 



Rf J-u R' Ml/ = d 



2[i/ 2 + f | +(d-l)H' 



R = d 



(d + l)H 2 + 2 



H 

N 



where we used 



a NH 2 H 

J-T _L 

N 2 a NN~ N' 



(4) 

(5) 
(6) 

(7) 



For simplicity we shall limit the discussion to a homogeneous scalar field 4>(t). Then 
the spatial volume can be integrated out from the measure in equation (JTJ), which we 
rewrite as 



S = J dtNa d [C g + Cf] 



The Hamiltonian constraint 5S/8 N\ 



N=l 



is 



' ' ON 

giving the Friedmann equation 

u\u — i 1 1 • rl ~ 

where 



C a + N d -±L - dHN d -h - N d dC * 



- P 



N=l 



dN dtdN 
d{d - l)FH' z — RF — f - 2dHF + 2(p + P4> + £p c 



\^ 2 + V, 



(8) 

(9) 
(10) 
(11) 



p c = 3a 4 4 - d[4 Cl EH 3 + (d - 3) Cl H 4 + c 2 (2EHH + 2dH 2 H - H 2 + 2HH)}. (12) 
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Here F = df/dR, E = and 

c\ = 2a% + da 2 + d(d + l)a 3 , (13a) 
c 2 = 4d! + (d + l)a 2 + 4da 3 . (136) 

Note that all three energy densities, p, p^ and p c , depend on the field <fi. In four 
dimensions (d = 3), the coefficients (jl3aj) - (J136j) read c\ = 2a\ + 3a 2 + 12a 3 and 
c 2 = 4(a! + a 2 + 3a 3 ), while the i7 4 term in p c vanishes. At low energy it was shown 
that the unique higher-order gravitational Lagrangian giving a theory without ghosts 
is the Gauss-Bonnet one (a x = a 3 = 1, a 2 = —4, a 4 = —1). In this case, c 2 vanishes 
identically while c\ = 2 + d(d — 3). With fixed dilaton coupling (5 = 0) equation (fT0|) 
reduces to the standard Friedmann equation in four dimensions, in agreement with the 
fact that the GB term is topological when d = 3. In three dimensions (d — 2), the GB 
higher-derivative contribution vanishes identically except for the 4 term. 

The continuity equation for the dark energy fluid contains a source term given by 
the coupling between this fluid and the string scalar field. We choose the covariant 



coupling considered in jlHl 113 HE] : SS p /5(p = —\/—gQ(4>)p, where S p = J d D Xy/—gC p 
and Q(4>) is an unknown function which we shall set to a constant later. In synchronous 
gauge we have 

p + dHp(l + w) = Q<fip, (14) 
while the equation of motion for the field <ft is 

lu(4> + dH<p) + V- gCf + 4a 4 ^ 2 (30 + dH<j) + S0) + ^0 - uo'^ - L ) = _Q P) (15) 
where the Lagrangian of the quantum correction is written as 



CP = d 



(d + l)ci# 4 + A Cl H 2 H + c 2 H 2 +a 4 4 . (16) 



Equations (JTUJ) — (JT3J) are the master equations of the physical system under study. 
2.3. Dilaton and modulus scenarios 

In order to reproduce general relativity at low energy, we impose that the leading term 
is given by the Einstein-Hilbert linear Lagrangian, i.e., f((j),R) = F((p)R. 

The massless dilaton field arises in the string loop expansion of the low-energy 
effective theory, its vacuum expectation value being directly related to the string 
coupling (e.g., 50 J. It is possible to consider the physical properties of the dilaton 
in two different frames (see, e.g., [JTJ ESI)- m the "string frame" (subscript S), the 
dilatonic action corresponds to 

F= - U j = e- 4 >, (17a) 
V = 0, (17b) 

£ = (17c) 
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where A =1/4, 1/8 for the bosonic and heterotic string, respectively, whereas A = in 
the Type II superstring. In fact, the above expression for £ is the tree-level term in the 
full contribution of n-loop corrections, given by 

n=0 

where Cq = 1. In this work we shall take only the tree- level term ()17c|) into account. 
Note that our function £(0) is dubbed — |c/A£(0) in |35] . 

One can transform to the "Einstein frame" (/ = R, subscript E) by making a 
conformal transformation 

9$ = FT£*gjS> = e#*gj£>, (19) 

so that Fe = oje = 1, while the term is modified accordingly (note that £e = is m 
four dimensions). Although the analysis can be performed in any frame, the equivalence 
principle is preserved only in the Einstein frame (in the string frame, the dilaton coupling 
to gravity results in a varying gravitational "constant"). One can shift from the string 
to the Einstein frame by noting that 

H E = e^2 I H s - — <— I , (20) 



20 

He = eo-2 



H s + 




D-2 D-2 \ D-2 



(21) 



where dots in the left (right) hand side denote derivatives with respect to the Einstein 
(string) time coordinate. Two frames coincide when the dilaton is fixed. 
In the string frame, the equations of motion read: 

d(d- 1)H 2 = (2dH -<j))<j) + 2e 4 'p + Xp c , (22) 
-Qe<V = + dH<j) + -4« + 2a 4 A0 2 (30 + dH<j) - 2 ) + ~ (0 2 + R^j , (23) 

where one has H = —0 in the definition of p c . 

In general, other moduli appear whenever a submanifold of the target spacetime 
is compactified with compactification radii described by the expectation values of the 
moduli themselves. In the case of a single modulus (one common characteristic length) 
and heterotic string (A = 1/8), the four-dimensional action corresponds to 

F = 1, (24a) 
u =3/2, (245) 
a 4 = 0, (24c) 

£ = - Al n [2eV(^)], (24a0 
Id 

where n is the Dedekind function and 5 is a constant proportional to the 4D 
trace anomaly and depending on the number of chiral, vector, and spin-3/2 massless 
supermultiplets of the N = 2 sector of the theory. In general it can be either positive 
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or negative, but it is positive for the theories in which not too many vector bosons 
are present. Again the scalar field corresponds to a flat direction in the space of 
nonequivalent vacua and V = 0. At large the last equation can be approximated 

as 



which we shall use instead of the exact expression. In fact it was shown in Ref. that 
this approximation gives results very close to those of the exact case. 

In any realistic situation both the dilaton and the modulus can appear in the 
action. However, we can consider the effect of each of them separately for the following 
reason. As already mentioned, a varying dilaton is strongly constrained by gravitational 
experiments. Although there is still a possibility to allow for a non-trivial dynamical 
contribution, one can assume to freeze it out at the minimum of its non-perturbative 
potential. Once the dilaton is fixed, gravity is not coupled to the surviving modulus field. 
Also, while in the modulus case one considers the heterotic string already compactified 
to d = 3, in the dilaton case we should start from the 26- or 10-dimensional action and 
then compactify down to four dimensions. To set d — 3 into the action is equivalent to 
consider a stabilized modulus (or moduli) while keeping the dilaton as a flat direction 
of the theory. 

3. de Sitter and inflationary solutions: fixed scalar field 

3.1. Preliminary remarks and geometrical inflation 

The search for de Sitter (dS) solutions of the higher-derivative equations of motion 
is important at least for two reasons. First, they can give useful insight as regards 
sensible implementations of string-motivated scenarios in the early-universe inflationary 
context [2U ES] • Secondly, stable dS solutions may avoid the infall of the late universe 
in eventual finite-time singularities. In order to take into account the case of inflation, 
we shall consider small but non-vanishing time derivatives of the Hubble parameter. We 
introduce the slow-roll (SR) parameter 

(27) 

while leaving the other SR parameters specified implicitly in the second-order relation 
H = 0(e 2 ). Since the evolution equation for e is at second-order in the SR parameters 
themselves, these are constant if small enough (slow- roll approximation). Pure dS 
solutions correspond to e = ■ ■ • = 0, where dots stand for other parameters of the 
SR tower. The Friedmann equation (jlUj) can be truncated at lowest order in the SR 
parameters by neglecting 0(e 2 ) contributions^ Let us study the case of a constant 
scalar field {(f) = 0) in order to present a preliminary example. 

§ Typical inflationary scenarios are achieved via the introduction of a dynamical scalar field. If we 
consider p in the action Q as the energy density of the inflaton, then it would be trivial to find quasi 



£ « £ cosh 0, 



(25) 




(26) 
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In general, (quasi) de Sitter solutions can exist for a non-linear f(R) as in the 
Starobinsky model [S3]. In the case of linear gravity, we can fix F — 1 = £ without 
loss of generality. If the scalar field vanishes or is constant with vanishing potential 
(p<p — = S), the Friedmann equation becomes a differential Riccati-type equation in 
H, 

4:C 2 dH + 2(d-3) Cl H 2 + (d-l) = 0, (28) 

where we have factored out the constraint giving the trivial Minkowski spacetime H = 0. 
Clearly there are no dS solutions for d = 1 and d = 3. At low energies the only invariant 
Lagrangian which is second order in the Riemann tensor and contains no ghosts is 
the Gauss-Bonnet invariant, which coincides with the Euler characteristic of the target 
manifold in four dimensions. As c 2 = 0, not even inflationary solutions would exist 
in the GB case. Nevertheless it is worth studying the case of c 2 7^ for the reasons 
explained in the introduction. Then the 4D inflationary solution is 

H(t) = ~, (29) 
0C2 

with C2 < 0. When 1 7^ d 7^ 3, there are both dS and inflationary solutions. The former 
is given by [22] H 2 = — (d — l)/[2(d — 3)ci]. If d > 3, then c\ must be negative, which 
excludes the GB case. The quasi dS solution for equation (|25|) with c 2 7^ is 



#(t) = J d } tanh 

W 2(3-d)£ Ci 



(d-l)(3-d)ci 



(30) 



8(i 2 £ocI 

which is obtained by fixing the integration constant so that H(0) = and by assuming 
Ci, C2 < (H(t) is positive and monotonic in this case). Changing the sign of c\ results 
in a periodic tangent behaviour. The Hubble rate (|3*U|) approaches a constant in the 
limit t — ► 00. To summarize, higher-derivative terms can induce an inflationary phase 
of pure geometrical nature. This unusual situation|| is anyway avoided in any dimension 
when considering a Gauss-Bonnet coupling. This might be another reason why to fix 
the coefficients q to their GB value. 

A stability analysis can be performed as follows. Let H c be a solution of equation 
(I28|) and H = H C + 5H be another solution given by the first plus a small time-dependent 
perturbation 5H 7^ 0. If the perturbation decays in time, then the solutions are stable. 
Dropping 0(5H 2 ) terms in equation (|25|) . one gets 

5H(t) = a(t) (3 - d)ci/(dC2) , (31) 

where we have used the definition of the Hubble parameter. Therefore the dS solution 
for d > 3 exists and is stable if ci,c 2 < O.^f Since the inflationary solutions can be 

de Sitter solutions. Here our aim is to consider the pure gravitational action with p = and see what 
kind of cosmological evolution is obtained in the presence of curvature corrections. 
|| Curvature-driven accelerating models were already considered in, e.g., |55|-[62| and references therein. 
See also for a case which can be compatible with solar system experiments. 

% This result, confirmed by the following numerical analysis, is in disagreement with [22], where it was 
claimed that dS solutions are unstable regardless of the parametrization. In particular, we were not 
able to reproduce their equations (33)-(36). 
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regarded as perturbations of the dS ones, intuitively these are stable when dS solutions 
are; indeed the inflationary solution (|3Uj) is stable. 

In the 4D case corresponding to equation f)29j) . the perturbation is constant and 
therefore the solution is not stable, as could be guessed by the fact that there is no dS 
background to perturb. 

We caution the reader that from the point of view of quantum gravity, models 
with a non-GB parametrization are unstable because of the presence of ghost fields. In 
this respect, whenever we mention the stability issue for such models we are implicitly 
referring to the classical stability against linear perturbations of a classical solution 
within a phenomenological effective theory. The details of the latter are still far for 
being established, as we shall discuss in the concluding section. 



3.2. Phase space analysis 

One may wonder if the slow-roll approximation we used in the previous subsection is 
really justified. A priori the H term may play a relevant role in the analysis, since the 
stability in that direction was not considered. In the next sections we clarify this point 
and generalize to the case of a dynamical scalar field. We shall consider the cases c 2 7^ 
and c 2 = (GB scenario) separately. 
By introducing the variables 

x = H , y = H , u = <f> , v = (j) , z = p , (32) 

we obtain the phase space evolution equations 

x = y, (33) 

(d-l)Fx F 2z + ujv 2 + 2V 3a 4 f 4 

y — 1 1 

4£c2 2£c 2 4d£c 2 x 2dc2X 

4ciHx 3 + (d- 3)cix 4 + c 2 (2Exy + 2dx 2 y - y 2 ) 

2c^c ' (34) 

ii — v , (35) 

-duxv - V + i'C^ - Aa^v 3 (dx + S) - cuv + (u'v 2 + F'R)/2 - Qz 
v = i 10 c 2 ' ( 36 ) 

z = [-dx(l + w) + Qv] z , (37) 
where 

£(?) = d[(d + l)cix 4 + Acix 2 y + c 2 y 2 } + a 4 t> 4 , (38) 

R = d[(d+ l)x 2 + 2y]. (39) 

In what follows we do not take into account the scalar field potential, which is not 
present in a perturbative string framework. In addition the coupling Q is assumed to 
be constant. 

In the GB case (c 2 = 0) differentiating the Friedmann equation ()10|) with respect 
to time gives the following equation 

2d[F +(d- l)Fx + \2dix 2 + 4(d - 3)c^x 3 ]x = (RF - /)" - 2dFx - d(d - l)Fx 2 
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+2z + uv 2 + 2uvv + 2vV + Qa^itv + 4fu) - 2d{d - 3)c^x 4 - 8dc^x 3 . (40) 

Let us search for a dS solution which is characterized by x = 0, y = with a 
vanishing fluid (i = z — 0). We find that 

2£d{d - 3)cix 4 + 8^CiSx 3 + d(d - l)Fx 2 + 2dFx - ujv 2 - 6a 4 £v 4 = . (41) 

When the scalar field (f) is fixed with vanishing potential (£ = 1 = F), equations (fTUjl 
and dlH) reduce to 

d(d - l)H 2 = 2p-2d (d - 3)dH 4 + c 2 (2dH 2 H - H 2 + 2HH) \ , (42) 
p + dH(l + w )p = 0. (43) 
Defining the vector X = (x, y, z) f , the above equations of motion can be written as 

x = y, (44) 
y = - dxy + a; H , (45) 

2C2 2x 4C2 2(3C2X 

i = — d(l + w)xz , (46) 

where we assumed that the perfect fluid has a constant barotropic index w. In the 
following we treat the case w = — 1 (cosmological constant) separately. Also, we consider 
a space with d > 3 in this section. 

3.3. Perfect fluid with w ^ -1 

The dS fixed point X c = for the previous set of equations is 



'• = ^W^' Sc = ' 2c = ' (47) 

where we have considered only the expanding case. When d > 3 we require the condition 
ci < for the existence of the dS solution. 

Linearizing the equations of motion under a perturbation 5X = X—X c with respect 
to a solution X c , one obtains 

5x = Sy, (48) 
3(3 - d)c 1 2 (d-1) z c 



8y 



2c 2 " Xc dVc 2x 2 r 4c 2 2dc 2 x2 



5x + I — — dxA 5y H fe, 



x c / 2dc2X c 

(49) 

5i = — c?(l + w)z c 5x — d(l + w)x c 5z. (50) 

This can be rewritten in compact notation as SX = M5X, where M is a 3 x 3 matrix 
with elements and eigenvalues 7 defined by the Jordan constraint 

7 3 - ("133 + m 22 )l 2 + (m 22 m 33 - m 2 i)7 + (m 21 m 33 - m 31 m 23 ) = . (51) 
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O.M> 




-0.60 I — ■ — ■ — ■ — ■ — i — ■ — ■ — ■ — ■ — i — ■ — ■ — ■ — ■ — i — ■ — ■ — ■ — ■ — i — ■ — ■ — ■ — ■ — I 
0.5 1 1.5 2 2.5 

X 



Figure 1. The phase space plot for d = 4, £ = 1, F = 1, w = —0.9, c\ = — 1 
and C2 = — 1 when the scalar field <\> is fixed. In this case the de Sitter fixed point 
(x c ,y c ,Zc) = (1.22,0,0) is a stable attractor. 



For the dS fixed point given by equation (jTfjl one has m 3 i = — d{l + w)z c = 0. In this 
case equation (foTj) is factored out as 7 = m 33 and 7 2 — m 2 27 — m%\ = 0. Therefore we 
find that the eigenvalues are 



71 = -d(l + w)x c , 72,3 = - 

This implies that the dS solution is unstable for a phantom fluid, since 71 is positive 
for w < — 1. When C2 < the real part of both 72 and 73 is negative. Therefore the dS 
expanding solution is stable for w > — 1 and Ci < 0, as is illustrated in figured When 
c 2 is positive we have 71 > and 72 < 0, implying that the dS solution is a saddle point. 
In summary, there is a dS solution for c\ < and d > 3 which is a stable attractor for 
C2 < and w > — 1. Note that parameters with this sign choice do not avoid either 
ghost-free graviton scattering amplitudes or naked singularity structures [HI]. 

3.4- Cosmological constant (w = —1) 

In this subsection we shall discuss the situation in which a cosmological constant A is 
present instead of a barotropic fluid. This corresponds to dropping equation and 
replacing p for A. Then we find that the dS solution satisfies 

A = d(d - 3) Cl x 4 c + ^_±K j yc = . (53) 

There exists one solution for this equation when (d — 3)ci is positive, whereas two 
solutions exist for (d — 3)ci < as long as the cosmological constant satisfies < A < 



-dx c ± \l d 2 x\ + 



2{d- 1) 



c 2 



(52) 
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d(d — 1) 2 /[16(3 — d)ci]. In the latter case each solution x = x c belongs to the range 



(d-1) , (d-1) 

0<x c <x M ^,/ i ^-^-, and *m<x c <^ w -^-. (54) 

One can compute the eigenvalues of the 2x2 matrix for perturbations Sx and Sy about 
the dS fixed points given by equation (JH3j) . They are 

71,2 = 2 (~ dx c ± ^Jd 2 x\ + 4m 21 ) , m 21 = [4(d - 3) Cl ^ + (d - 1)] . (55) 

If (d — 3)ci is positive, we find m 2 i < for c 2 > 0, which means that both 71 and 7 2 are 
negative. Then this case corresponds to a stable attractor. Meanwhile the dS solution 
is unstable for [d — 3)ci > and c 2 < 0, since one has 71 > and 7 2 < 0. 

When id — 3)c± < there exist two dS fixed points given by equation (|54J) . For the 
first critical point in equation (|54*|) . one has m 21 < for c 2 > and m 2 i > otherwise. 
Therefore the dS solution is stable for c 2 > and a saddle point otherwise. Similarly 
we find that the second critical point in equation (|H4*jl is stable for c 2 < and a saddle 
point otherwise. 

From the above argument it is clear that the stability of the dS solutions crucially 
depends upon the signs of C\ and c 2 . 

4. Linear dilaton case 

From now on we shall study the cosmological evolution in the presence of a dynamically 
varying scalar field. For the dilatonic action in the string frame, the condition (|4~T|) for 
the existence of the dS solution yields 

Ad(d - 3)cix 4 - 4d\civx 3 + d(d - l)x 2 - Idvx + v 2 - 3a 4 Aw 4 = . (56) 

This suggests that dS solutions may exist when v = <fi is a non-zero constant. 
Using v = in equation (jSBj), one gets 

— \d(d + l)cix 4 — AdXa^x — d(d + l)x 2 + 2dt> x — v 2 + Sa^Xv 4 = . (57) 

From the above two equations we find that the dS solution satisfies 

2Xdx 3 + 2X(c 1 vx 2 + a 4 v 3 ) +x = 0. (58) 

We note that the Minkowski solution (x = and v = 0) also satisfies equations (|56|) 
and J5ZJ). 

When d = 3 equations and (|55jl give 

_ -^ 2 (9Aa4^ 2 + 1) , , 

~ 6(l + 2A Cl t; 2 ) • 1 J 

For the bosonic correction with a GB term one has A = 1/4, 0,4 = —1 and c\ = 2. 
Then by combining equations (|58j) and (fo^j). we obtain the following dS solution with a 
linearly changing dilaton: 

x c = 0.62, u c = 1.40. (60) 
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Figure 2. The phase portrait (plot of y versus x) of the system in the linear dilaton 
case with d = 3, c\ = 2, A — 1/4, = —1 and C2 = — 1. Trajectories starting anywhere 
in the phase space converge at (x c , y c , z c ) — (0.62, 0, 0) for a constant value of v = 1.40. 
The fixed point corresponds to a stable de Sitter solution. 

This agrees with the result in Ref . [H] . It is clear from equation (JoTJj) that the dS solution 
does not exist for a 4 = for the GB coupling (ci = 2). If we allow for a negative value 
of c\, it is possible to have a dS solution even when 04 = 0. Note that these results do 
not depend on c 2 . 

The issue of the stability of such solutions is presented in more details in the 
Appendix (including the case of C2 = 0). A rather counter-intuitive result is the 
following. Although the dS solution is insensitive to the actual value of c 2 , its stability 
depends upon the sign of c 2 . More precisely, the solution is stable against linear 
perturbations for 04 < 0, c\ > 0, c 2 > 0, whereas for negative c 2 one of the eigenvalues 
has a positive real part and the solution is unstable. We have checked other situations 
for d = 3, which can be summarized as follows (we considered only real solutions with 
x c > 0): 

• Ci > 0, a 4 > 0: There is no real solution. 

• Ci > 0, a 4 < 0: Real solutions (a 4 < — (9XV 2 )" 1 ) are stable only for c 2 > 0. 

• Ci < 0, a 4 > 0: Real solutions (ci < — (2Af 2 ) -1 ) are stable only for c 2 > 0. 

• ci < 0, 04 < 0: Real solutions are stable only for c 2 < 0. See figure El 

We should stress that the dS solutions we found appear in the presence of the in- 
dependent terms. In the low-energy regime where the effect of higher curvature terms 
is neglected, we obtain a logarithmic evolution of the dilaton corresponding to the 
pre/post-big-bang solution in pre-big-bang cosmology (121113] ■ Later we shall briefly 
discuss the case of a logarithmic dilaton in the string frame. 
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4-1. Impact of the barotropic fluid 

In the case of vanishing coupling Q, the above stability analysis holds as long as we 
assume a standard equation of state characterized by w > — 1. In the presence of a 
phantom fluid (diverging z c ) the eigenvalue 72 = —d(l + w)x c is always positive for 
x c > and all solutions are unstable (see the Appendix). If Q 7^ the dS solutions are 
stable when 



V c 

where we assumed v c > 0. Therefore the phantom fluid with positive Q is always 
unstable. 

4-2. Cosmological constant (w = —1) 

When w — — 1 and z = A = const, the left-hand side of equation (JHH)) has an extra term 
— 2Ae u . This term vanishes or diverges depending on the sign of u. In order to find a 
non-trivial dS solution with cosmological constant, one has to assume that the dilaton 
field u approaches asymptotically, v — ►(), u — > as t — > 00. In this regime, equation 
(IHTjl gives 



which is positive for ci < 0. The value of the cosmological constant is fixed by the choice 
of the coefficients. We ran our numerical code for various values of c\ < and C2 and 
checked that the dS solution in four dimensions is unstable against linear perturbations. 
Note that this behaviour does not coincide with what has been found for a fixed dilaton 
in section 13.41 This is because the dynamical equations are affected by the presence of 
the scalar field even before reaching the asymptotical regime. 

A quick inspection of the dynamical equations shows that in the Einstein frame 
there are no dS solutions, except under the same assumption adopted for the 
cosmological constant case (u —>■ asymptotically). 

5. Logarithmic scalar field and the late time cosmological evolution 

In the first part of this section we shall consider the cosmological dynamics in the 
presence of the compactification modulus under the assumption that the dilaton is 
fixed. This means that the analysis is identical both in the Einstein and string frame. 
In section l5~5l the logarithmic dilaton in the string frame will be studied as well. In order 
to keep our discussion as general as possible, we shall not fix the dimension d of the 
target space except at the end, where the four-dimensional case d = 3 will be inspected 
among the others. The equations of motion for the modulus action corresponding to 
Eqs. fl2H-(|243) read 




(61) 



x = y, 



(63) 
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2z + 3v 2 /2 - d(d - l)x 2 - 8d Cl £x 3 - 2d(d - 3)c^x 4 - 2dc 2 £y{2Zx + 2dx 2 - y) 



y 



y£'[(d + l) Cl x 4 + A Cl x 2 y + c 2 y 2 ] 



Adc 2 ^x 



dxv z, 

3 ' 



(64) 
(65) 
(66) 
(GB 



z = [—dx(l + w) + Qv] z. 

While only derivatives of £ appear in the equations of motion for d = 3 and c 2 
case), there are non- vanishing contributions of £ itself for general coefficients q. When 
c 2 = 0, the equations of motion for x and v read, from equations (jjQj) and (}36|) . 

z + 3w/2 - dcix 3 [4£ + (d - 3)ix] 



x 



dx[{d - 1) + \2dix + 4(d - 3)c^x 2 } ' 



4 = 2 4^-^ + l)x 2 + 4i] - d X - - 



(67) 
(68) 



d(d - 1)±- - ^ - 2- + 2d Cl ^[4S + (d - 3)x] 



0. 



(69) 



while the Friedmann equation is 

x 2 3 „ z , , 
t>^ 2 f * " v z 

When finding solutions one can use the equations for c 2 ^ and then set C2 = at the 
end, although equations (fB7 j) - (fBnj) are required for numerical purposes. In addition c\ 
can be set equal to 1 and absorbed in the definition of £o, so that the coefficient c 2 is 
the only free parameter of the higher-order Lagrangian. 



Following the approach of 



we search for future asymptotic solution of the 



form 



x 



u 



uit?, y 
u + lu 2 lnt, 



z ~ z t QuJ2 exp 



z t c 



0J 2 

V ~ — . 

t 



d(l+w)u! 0+1 
(3 + 1 

-- -1, 



where the barotropic index w is constant and 
a = Qu 2 — d(l + w)uj\. 

We define 



(70 a) 
(706) 

(70c) 

(70 d) 
(70 e) 

(71) 

(72) 



so that in any claim involving the sign of 5 (5) a positive C\ coefficient is understood. 
In order to find a solution in the limit t — > +oo, one has to match the exponents of t to 
get algebraic equations in the parameters j3,u>i,c 2 and 5. 
We shall consider the following four cases: 



(i) A low-curvature regime in which £ terms are subdominant at late times. 
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(ii) An intermediate regime where some terms in the equations of motion, either coupled 
to £ or not, are damped. 

(iii) A high-curvature regime in which £ terms dominate. 

(iv) A solution of the form (j70a|) - (j70e|) for the full equations of motion. + 

Our results will be in agreement with [22] (and [21], taking into account our extra 
a/2/3 factor in the normalization of v ) for d = 3 and z = 0. 

In order to keep contact with observations, we shall compare the obtained values 
of tO\ of the geometrical (z = 0) modulus solutions with the estimate of the deceleration 
parameter from supernovae data jS], which does not depend on the Einstein equations 
governing the background evolution [HE] (while the relation between the deceleration 
parameter and the dark energy barotropic index does). Solutions in agreement with 
observations might be used to describe the recent evolution of the universe. We shall 
not consider the logarithmic dilaton solution since it does not belong to the runaway 
dilaton models of fTJ ^] and may be in conflict with gravity experiments, although we 
have not checked this explicitly. 

As we will see later, one has /3 — — 1 for all except one case. Then u\ is given in 
terms of the deceleration parameter q = e — 1: 

<* = IT? < 73 > 

hat onnafinn \ nil r»an rlocr*riK^ nnKr nncmnlnrripc ^snt h fi — 1 ^ 



The first thing to note is that equation (|70a|) can describe only cosmologies with q > — 1 
that is, the effective phantom regime is avoided. 

At late times (redshift z ~ 0), the observational data give the constraint — 1.0 < 
<?o ^ —0.5 6J, where q = q(z)\ z ~o. This results in the condition for LO\ 

> 2, (74) 

if the obtained solution describes the present universe. 

Of course this constraint is loosened if we use the 3a observational bounds, but we 
need to exclude the region q < — 1 for the above argument. For —1.3 < q < —0.2 (after 
exclusion of the extremal negative region), one has ^1-2. 

In order to allow a super-accelerating phase (q < —1), we can shift the origin of 
time so that t — > t — t s . For positive increasing time t <t s , the solution 

x = (75) 

t-t; v ; 

+ Low- and high-curvature regimes may be viewed even as "weak coupling" and "strong coupling" 
models, respectively. These correspond to (late-time) stages at which higher-order terms become either 
subdominant or dominant because of the relative magnitude of the theory's coefficients, evaluated at 
some moment of the cosmological evolution. Quotation marks are in order since actually the curvature 
is not related to the modulus expectation value u; however, in the Friedmann and modulus equations the 
higher-order term £p c dominates both when £ 3> 1 and p c S> 1, and the two descriptions are equivalent 
(with a 4 = 0). The varying coupling picture may be confusing in our context, since curvature terms 
may dominate because of their asymptotical time dependence, for natural values of their couplings, 
while the weak/strong coupling approximation is assumed to hold in a time interval At where the 
dynamical variables have the above time dependence. 
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is expanding when oj\ < 0. The constraints on the parameters can be found by 
considering the previous asymptotic behaviour for t — > 0, that is near the singularity at 
t s . Then the previously discarded 3a region —1.3 < go ^ — 1 corresponds to 

< -57.8. (76) 

In order to satisfy this bound, however, the parameters of the solutions will require 
rather unnatural adjustments. 

We shall not discuss the stability of these solutions since they will hold at times 
t ~ t s , that is, only near the sudden future singularity. 

5.1. Structure of the solutions and behaviour of the perfect fluid 

Before proceeding, let us explain the general structure of the solutions. Each of the 
equations of motion, under the substitution (j7()qj) - (j7()e|) . is of the form 

E(^+^ r ) =0, (77) 

i 

where gf = gf((3, u>i, uj 2 , C\, c 2 , £o> • • •) are the functions of the parameters of the theory, 
and af = at(P, tui, u) 2 , Q, w) are non-negative exponents. In the limit t — > oo, an 
asymptotic solution with 

a+ = 0, a~ > 0, (78) 

is Yliidt = 0- If the gt coefficients do not depend on £ , then the solution represents 
a low-curvature regime. If dg^ /d£o ^ for some j, then we will say that the solution 
describes an intermediate regime, when some (but not all the) higher-order terms 
survive. If dgt /d£,o ^ for all i, only the higher-order terms contribute and the solution 
is in a high- curvature regime. 

In the limit t — > 0, the solution (fT8|) becomes J2i9i^ = 0> with eventually (in fact, 
only in the intermediate regime) some subdominant g~t~ a * terms dying away. Therefore 
all the above cases are interpreted in the complementary way, and low (high) curvature 
solutions become high (low) curvature. 

The solution satisfying at = a^ = is valid at all times and the only distinction 
between solutions with or without future singularity is the sign of oj\. 

As regards the contribution of the barotropic fluid, when (3 = —1 we can 
parametrize its time dependence in equation (|77jl as 

g z (t) = z t a - a *, (79) 

where a is defined in equation (fTTj) and a* is a constant determined by the kind of 
solution (low/high curvature, exact, etc.) one is searching for. 

We can adopt two perspectives. In the first one, we do not fine-tune the fluid 
parameters and impose that g z (t) vanishes asymptotically (purely geometrical solution). 
In this case, z — ► asymptotically in the equations of motion and a* — a G {ctj 7 }- From 
equation (|75jl it follows that a < a*. In the second one, the expansion is driven by the 
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perfect fluid and g z = const, that is, a — a* G {af }, and a = a*. We summarize these 
constraints in the damping/matching condition 

a < a*. (80) 

When the inequality holds, the fluid decays away in the limit t — > +00 and is consistently 
discarded. 

When equation (fKUJ) is an equality, the fluid contributes to the asymptotical 
equations together with other geometrical terms and affects the dynamics of the system. 
However, this does not mean that the energy density z(t) increases in time, since a can 
be negative. In this case the barotropic fluid can be regarded as a model for dark energy 
with 

1 d(l + w eS ) 



ux 2 



1 



_ 2w - Qu 2 

w °« = -JTcm- (82) 

Then constraints on q would result in bounds for w e g. In the asymptotic limit t — > +00, 
solutions exist (and do not burst out) only for the effective barotropic index w e g > — 1. 
However, in this section we shall constrain only the geometrical solutions with a scalar 
field 4>, postponing a more complete discussion on dark energy to the next section. 

To keep the notation clear, we shall write explicitly the time dependence of fluid 
terms and discuss separately the damping and matching cases. In the last section we 
shall compare numerically evolved dark energy scenarios to the analytical solutions 
below. In some cases the fluid will be dominating or decaying and either the matching 
or the damping condition will hold at the infinite future, with fixed constant parameters 
uj\ and UO2. At intermediate times, where avoidance or incoming of a big rip become 
manifest, there is enough freedom in varying w and Q to span a wide range of dark 
energy phenomena, since dg z /dt 7^ in general. 

5.2. Low- curvature and intermediate regimes 

The solution in a low-curvature regime corresponds to the one where ^-dependent terms 
are neglected. The asymptotic future solution we consider is (d ^ 1) 

P = -1, u 2 < 2, (83) 

together with the constraints 

1 2Qz r+ 2 

U1 = 1 — ^ — ' ( 84 ) 

a Sdu>2 

3w| = 2d(d - - Az t a+2 , (85) 
a < -2, (86) 

so that the energy density z(t) decreases in time. £ terms decay away regardless of the 
value of c 2 . According to equation (|S5jl. if oj 2 > and the barotropic fluid is phantom- 
like (w < —1), then Q should be negative and of high enough magnitude in an expanding 
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universe (oji > 0).* The reader can check other situations. In the damping case a < —2, 
the solution is U\ = 1/d and cu 2 = ±^/2(d - l)/(3d). This case (for d — 3) was labelled 
in Ref. [23! • 

In the marginal case a = —2, the scale factor expands for Q/uo 2 < 3/(2zq). This 
requires either an ordinary fluid with Qu 2 > —2 or a phantom fluid with Qu 2 < —2. 
When Q = one has w — 1. In addition equation (JHo")) states that uj\ > 2zo/[d(d — 1)]. 
We have numerically checked that this solution exists for suitable values of Q. 

Another solution, valid in an intermediate regime and only for the c 2 = and 
uoi > case, is 

P = -2, u 2 = 5, Q < -2/5, u\ = — ^ (15 - 2Qz t 5Q+2 ) , (87) 

16d5 

for a non- vanishing fluid. We require the condition 5 > in order to have an expanding 
scale factor a(t) ~ a exp(— Ui/t). Actually this goes from to a , that is, it reaches 
Minkowski spacetime asymptotically. If the fluid decays, then one recovers the 
solution of Ref. j2H] with d = 3 and z = 0. Interestingly, the coupling Q cannot vanish 
when z ^ 0. 

For both the low-curvature and intermediate regimes, there exist no solutions 
corresponding to (j70aj) - (j70ej) in the cosmological constant case (w = —1), since the 
A/v 2 term diverges asymptotically. 

We performed numerical simulations for c 2 = 0, z = 0, and u 2 > with general 
initial conditions determined by the Hamiltonian constraint. With either positive 
or negative values of 5, we found that in the asymptotic future the solutions tend 
to approach the general-relativistic (GR) solution rather than the intermediate 
solution (JBZJ). One example is shown in figure El corresponding to the case without 
barotropic fluid and d = 3. The asymptotic solution is actually described by constant 
values of u)\ and uj 2 characterized by equations an d (jHHjl . i.e., w\ = 1/3 and 
uj 2 = 2/3. The 5 < case is consistent with the result in Ref. [23], where non-singular 
cosmological solutions were constructed by joining the future GR solution (J83|) with the 
past intermediate solution (|S7jl . In fact we numerically found that the asymptotic past 
solutions are well described by equation ()87|) for negative 5. One can either choose a 
positive or negative sign of </> in order to obtain non-singular cosmological solutions. 
When 6 > 0, we found that the past solutions are not described by equation (|57jl. which 
is associated with the fact that singularity- free solutions do not exist in this case 

It is easy to check the stability of low-curvature solutions without fluid analytically. 
Using equations (|57j) -(|5U |) for z = 0, we obtain, after dropping ^-dependent terms, 

x w -dx 2 , (88) 

which gives 

d* = (89) 

* Constraints on the sign of Q will depend on the sign of ui2- Assuming v < would simply flip sgn(Q) 
in these considerations. 
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Figure 3. Evolution of lo\ and lo 2 for c 2 — 0, d = 3, z = 0, and £o = 2 without a 
barotropic fluid. We find that the system approaches the low-curvature solution 1831) 
characterized by constant lo\ and u> 2 {u>i = 1/3 and lo 2 — 2/3). 



Considering perturbations 8u\ (t) around the critical point ui^ = 1/d, we obtain 



(90) 



Since v and x are proportional to each other, the perturbations in u> 2 decay as 8u>x, 
thereby showing the stability of the GR solution. In the presence of a fluid with Q = 0, 
it can be shown that the solution is stable again, provided that w > —1. 
The solution with a sudden future singularity at t s corresponds to 

(3 = -I, u 2 >2; (91) 

and an eventual cosmological constant term damps away. The parameter constraints 
are equations (fMj) and (f£K|h In the marginal case 

a = -2, (92) 

the scale factor expands for Q > 3u2/(2z ) > 3/zq. This requires a phantom fluid with 
w + 1 = (2 + Q0J2) / {duj\) < 0. We have numerically checked that this solution exists for 
suitable values of Q. 



5.3. High curvature regime 

When higher-order terms become important we can consider the asymptotic future 
solution 



/? = -!, 



U)2 > 2, 



a < u>2 — 4, 



(93) 



Dark energy and cosmological solutions in second-order string gravity 



22 



together with constraint equations 

5dul[(d + l)ul - 4w x + c 2 ] - Q^™ +4 = 0, (94) 
d5cuf[(3 - d)uj\ - Auj 2 uJi + c 2 {2du 1 + 2cu 2 - 3)] + z t a ' W2+A = 0. (95) 

In the GB case (c 2 = 0) the solution corresponding to a decaying fluid (a < u 2 — 4) is 

4 3-d . , 

Ul = d + 1' " 2 = d + 1' (%) 
which contradicts the condition (|93|) in any dimension. Therefore in the GB scenario 
with ui 2 7^ 4 only the marginal case a = uj 2 — 4 is allowed. 
When 

uj 2 = 4, (97) 

the cosmological constant term (w = —1) survives; an additional fluid term would be 
trivial since a = 0. The constraints of this solution are 

(d+ 1)^-4^ + 02 = 0, (98) 
d8ul[{3 - d)uj\ - 16wi + c 2 {2du 1 + 5)] + A = 0. (99) 

In the GB case (c 2 = 0) we find 

^ = JT~v 4 (100) 

A =5d(5d+l) ■ (101) 

\a+ I J 

The presence of a non-vanishing cosmological constant is mandatory for the consistency 
of this solution. For example one has A = 485 for d = 3. It is remarkable that equation 
(|97|) is the only situation where a non-zero cosmological constant is not only allowed 
but even required in the GB case. 

In order to check whether the solutions really approach the one given by 
equation ([9*7|). we run our numerical code for d — 3, c 2 = 0, and 5 > with A = 485. 
As illustrated in figure IH the future asymptotic solutions are described by equations 
(|97|) and ()100|) with constant u)\ (=1) and u 2 (=4). We chose several different initial 
conditions and found that the solution (|9*7j) is a stable attractor. We also considered the 
case of a negative cosmological constant, but the solution (|9*7j) is found to be unstable. 

When c 2 ^ 0, not all values and signs of the coefficient will give real expanding 
solutions. Since 



^ - 2± V 4 -( rf + 1 ) C2 (102) 

the reality condition reads c 2 < A/ (d+1). If c 2 < 0, only the solution is expanding, 
whereas both solutions are for < c 2 < Ajid + 1). Again, the sign of the cosmological 
constant will depend upon 5. 

Numerically, the solution of this system for d = 3 approaches the analytical value 
of uo^ (and the associated uj 2 ) for positive A. The u± branch is found to be unstable 
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Figure 4. Variation of w\ and U02 for C2 = 0, d = 3, and £0 = 2 in the presence of 
a cosmological constant given by A = 48<5. The asymptotic behavior of H and 4> is 
characterized by p = — 1, uj\ = 1 and LU2 — 4 as estimated analytically. 



regardless of the sign of the cosmological constant. Therefore the only stable solution 
is oof with A > 0. 

There are non-trivial positive solutions even with A = 0, corresponding to the 
critical case 

2 4 



Ul = dTT' C2 = J+i' (103) 

We chose several different initial conditions and found no stable attractor. Since we did 
not check it for an extensive set of initial conditions, we cannot claim anything definite 
regarding the stability of this solution. 

The sudden future solution corresponds to 

[3 = -1, u 2 < 2, (104) 

which describes a high-curvature regime where the ^-dependent terms are dominant 
(a cosmological constant term is suppressed). The constraints on the parameters are 
equations (jHlJ), (JUSl) and 

a > uj 2 - 4, (105) 

where z — > when the last equation is an inequality. In this case the solution is 
contracting for c 2 = 0. Conversely, for z = and c 2 7^ the only expanding solution is 
u>i given in equation ()102j) for c 2 < 0, together with 

uj 2 = 3 — du^ . (106) 
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5-4- Exact solution 

An exact solution which is valid at all times is 

= -1, cu 2 = 2, (107) 
together with the constraints on w\\ 

25diu 2 [(d + l)cuf - 4wi + c 2 ] - §du x + 6 - Qz t a+2 = 0, (108) 
25duf[(d - 3)w? + 8ut - c 2 (2duj x + 1)] + d(d - l)u}\ - 6 - 2z t a+2 = 0, (109) 

and equation (jHEJ). The cosmological constant term is forbidden in this case. 

Let us consider the case with a vanishing fluid (z — > 0). After some manipulations, 
the above equations can be written as 

A5ujl\(d - + 2 - dc 2 ] + {d- l)wi -6 = 0, (110) 
25cji[d(rf+ l)^ 2 - 2(d + + 4 - dc 2 ] - (5d + 1) = 0. (Ill) 

When c 2 = 0, these equations completely fix the parameters of the theory. In four 
dimensions, for instance, the only real root of the above equations is 

{ut, 5) = (0.21,13.53), (112) 

in agreement with Ref. (solution B^). Note that the original coupling 5 is not 
fine-tuned, although 6 is, since we have a freedom in choosing uq. There is no solution 
with a sudden future singularity. 

Numerically we found that the exact solution (|107|) is not stable. As we already 
pointed out in section 15. 2\ the future stable attractor for A = corresponds to the GR 
solution given by equation (JH3~j) for both positive and negative 5. The asymptotic past 
solution is described by equation (j%7|) . with 5 < 0. 

In the case c 2 ^ 0, from equation ()110|) one can show that a solution of the system 
is given by (d ^ 1) 

6 8 ~ (d-lf , , 

We checked numerically that this solution is unstable. 

Other (real) solutions, with either uo\ > or io\ < 0, can be found by fixing the 
parameters known a priori from the theory, that is, c 2 and 6. 



5.5. Logarithmic dilaton 

With little effort, we can check whether equations (j70aj) - (j70ej) are a solution of the 
dilatonic model in string frame for t — > +oo. It is easy to find that the only possible 
case is (3 = —1, corresponding to a low-curvature regime where ^-dependent terms are 
damped away. Then the model reduces to a scalar-tensor theory in Einstein gravity. 

In the presence of a barotropic fluid, the matching/damping condition reads 
a < —2 — u 2 , i.e., 

(l + Q)oo 2 -d(l + w)uj 1 + 2<0, (114) 



Dark energy and cosmological solutions in second-order string gravity 



25 



while the constraints on the parameters read, with the above matching condition, 

du\ - dw x + u 2 + (1 + Q)e u °z t a+UJ2+2 = 0, 

d{d - - 2du 1 u) 2 +oj 2 2 - 2e u °z t a+UJ2+2 = 0. 
When a = — 2 — u 2 , the solutions are 



±d Yd 



1 - 2(1 + Q)e"°z + v/l+4(l-Q)e«o Zo 
^{-1-^1 + 4(1-3)^0 



(115) 
(116) 

(117) 



± k 2d 



1 - 2(1 + Q)e u °z + y/1 + 4(1 - Q)e u °z 



~.± 



UJ 



l 



1 - 2(1 + Q)e u °z - VI + 4(1 - Q)e u °z 



'1181 



;ii9) 



1 +4(1 -Q)e u °z 



± \ 2d 



1 - 2(1 + Q)e u °z - v 7 ! + 4(1 - Q)e"°z 



Without fluid (Q = = z ), the expanding solution reads 

1 

v 7 ^- 1, 



(120) 

(121) 
(122) 



in agreement with Ref. [52]. One has ^0.58 and o; + ~ 0.73 for d = 3. This solution 
is found to be stable numerically. 

In the presence of a cosmological constant [zot a = A and Q = in equations ()115|) 
and ((HE)) ], we get uj\ = —2/(l + d) and c<j 2 = —2, which describes a contracting universe. 

The solution with a sudden future singularity at t s and (3 = — 1 corresponds 
to a high-curvature regime where ^-dependent terms are dominant. In this case the 
constraints read 

3Aa 4 c4 - \du>l[(d - 3)ciw? - Acmj^ + c 2 (2uj 2 - 2cLji + 3)] + 2e u ° z t a+UJ2+4 = 0, (123) 

2Xa^(du 1 - w 2 - 3) + - {dw?[(d + ^i^i ~ 4 ci^i + c 2 ] + a 4 w|} + Qe w °2 t a+a;2+4 = 0. 

(124) 

When A ^ z 7^ 0, the damping/matching condition for the fluid is 

(l + Q)cj 2 -d(l + w)cj 1 + 4> 0. (125) 

When the fluid vanishes or decays, we verified that there exist non-trivial solutions for 
suitable choices of parameters, both without and with a cosmological constant (z — A 
and u 2 = —4 in the above equations). For instance, in four dimensions and in the GB 
case {d — 3, c\ — 2, c 2 = 0, a 4 = —1), the only expanding solution is uj\ = —2.87 and 
A = —188.88. One can find other solutions by varying the parameters ci, c 2 , a 4 , A. 
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5. 6. Constraints from the current universe 

We now compare the observational constraints on uj\ for the recent evolution of the 
universe with the modulus solutions found in the previous section (d = 3). We consider 
the situation in which a perfect fluid is vanishing asymptotically. The results are 
summarized in Tabled with the 68% CL bound. We discarded the intermediate regime 
or Minkowski solution. 



Solutions t — > oo 



CO 



(0) 



(68% CL) Stability 



Low curvature 

High curvature, C2 = 

High curvature, c 2 ^ 0, A = 

High curvature, c 2 / 0, A / 

Exact, C2 = 

Exact, C2 7^ 0, equation ()113|) 




Yes 
A > 

uf, A > 



Table 1. Constraints on modulus solutions in the asymptotic future for the Hubble 
parameter H = cjit -1 , vanishing fluid, and d = 3. Blank entries are excluded by 
experiments or numerical analysis. 



We find that the logarithmic modulus solution with ghost-free GB parametrization 
and no extra fluid does not provide a viable cosmological evolution in the current late- 
time universe. However, we will see in the next section that the GB case in the presence 
of a cosmological fluid may display interesting features as regards the future evolution 
of the universe. At the 99% confidence level, the low redshift constraint on C2 for the 
high curvature, A ^ solution can be relaxed up to c 2 < — 1. Therefore we have shown 
that there are models which can in principle explain the present acceleration without 
using the dark energy fluid. 

The situation becomes more complicated in the presence of a barotropic fluid. 
The low-curvature solution can describe the very recent universe if Q is negative and 
non-vanishing. The other cases crucially depend upon the interplay between all the 
theoretical parameters. 

As regards the models with a sudden future singularity, the bound equation (fTE|) 
generally requires severe fine tunings for the parameters of the solutions. For instance 
one needs the condition Qzq > 10 3 for the low-curvature solution with u> 2 = 4 and d = 3. 
In the high-curvature case (|104|) with a vanishing fluid, we require that \c 2 \ > 10 4 . In 
the exact case (c2 ^ 0), \8\ < 1CT 6 . In the dilaton case with a cosmological constant, 
|A| > 10 7 - 10 9 . 

However, we found a remarkable theoretical result. This kind of models with c 2 = 
(Gauss-Bonnet case) always requires a non- vanishing contribution from the barotropic 
fluid for consistency. In the low-curvature case, even a positive coupling Q is necessary. 
In the dilaton case, a non- vanishing cosmological constant must be present if z = 0. 
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6. Gauss— Bonnet modulus gravity and the dark energy universe 

In the previous section we compared cosmological solutions obtained in second-order 
gravity with the present evolution of the universe. In this section we shall adopt another 
perspective, namely, we will investigate the phenomenology of the full equations of 
motion as models of future evolution of the dark energy universe. At the infinite future, 
the properties of such models will agree with some of the above solutions. 

In the typical phantom scenario (w < —1), the energy density of the fluid continues 
to grow and the Hubble rate eventually exhibits a divergence at finite time (big rip). 

It is expected that the effect of higher-order curvature corrections becomes 
important when the energy density grows up to the Planck scale. In fact, it was shown 
in |T3J inni EIZ| that inclusion of quantum curvature corrections coming from conformal 
anomaly can moderate the future singularities. 

Here we would like to consider the effect of O(a') string quantum corrections when 
the curvature of the universe increases in the presence of a phantom fluid. We shall 
concentrate on the modulus case with £ given by equations (|23|) and (|26|). Our main 
interest is the cosmological evolution in four dimensions (d = 3) in the presence of a 
GB term (ci = 2, C2 = 0) with a 4 = 0. We assume that the dilaton is stabilized, so that 
there are no long-range forces to take into account except gravity. 

From the discussion in section we found that the growth of the barotropic fluid 
is weaker than that of the Hubble rate when equation (|86|) is satisfied as an inequality. 
This condition is not achieved for a phantom fluid when the coupling Q between the 
fluid and the field is absent (Q — 0). 

We ran our numerical code by varying initial conditions of H, cf) and p in the full 
equations of motion. When 5 < 0, we numerically found that the solutions approach a 
big rip singularity for Q = and w < — 1 (see figure EJ). Meanwhile the condition (|86|) 
can be satisfied for negative Q provided that u 2 is positive. In figure El we show the 
evolution of H and p for Q = —5 and w — —1.1. In this case p decreases faster than t~ 2 , 
which means that the energy density of the fluid eventually becomes negligible relative 
to that of the modulus. Therefore the universe approaches the low-curvature solution 
given by equation (J55|) at late times, thereby showing the avoidance of big rip singularity 
even for w < —1. By substituting the asymptotic values u\ = 1/3 and 0J2 = 2/3 in 
equation (|5Sj). the condition for decaying fluid reads Q < 3(w — l)/2 = —3.15. We 
checked that the big rip singularity can be avoided in a wide range of the parameter 
space for negative Q. These results do not change even for smaller values of S such as 
\8\ = 0(1) (corresponding to |£o| = 0.1). 

When Q is positive, the condition (|KH|) is not fulfilled for 0J2 > 0. However our 
numerical calculations show that <fi becomes negative even if <fi > initially. We 
found that the system approaches the low-energy regime characterised by uj\ = 1/3 
and oj 2 = —2/3. Since co> 2 < 0, the big rip singularity can be avoided even for positive 
Q. In fact we numerically checked that the Hubble rate continues to decrease as long 
as the condition (jHU|) is satisfied in the asymptotic regime. 
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Figure 5. Evolution of H and p with £o = —2, w = —1.1 for (a) Q = and (b) 
Q = —5. We choose initial conditions as Hi = 0.2, <\>i = 2.0 and pi = 0.1. The curves 
(al) and (bl) represent the evolution of H for Q = and Q = —5, respectively, while 
the curves (a2) and (b2) show the evolution of p for corresponding Q. We find that 
the big rip singularity is avoided when the negative coupling Q is present. 

When 5 > 0, there is another interesting situation in which the Hubble rate 
decreases in spite of the increase of the energy density of the fluid. This corresponds 
to the solution in the high- curvature regime in which the growing energy density p can 
balance with the GB term (p m 24H 3 !; in the Friedmann equation). One example is 
illustrated in figure the big rip does not appear even when w < — 1 and Q = 0. 
Thus the GB corrections coupled with a scalar field <fi provides us several interesting 
possibilities to get rid of this singularity. 

By examining the condition (|80p. one may think that the big rip singularity may 
appear even for w > — 1 in the presence of the coupling Q. However we found that 
this is not the case. We show one example in figure corresponding to w = —0.8 and 
Q = +5. It turns out that <ft becomes negative even when <p > initially, as seen in 
the figure. We also checked that this property holds even for negative Q. This suggests 
that the system tends to evolve so that the condition 

a < when Q + 0, (126) 

is satisfied. Then the presence of the coupling Q does not lead to the big rip singularity 
for a non-phantom fluid (w > —1). 

7. Conclusions 

In this paper we have studied cosmological solutions in the presence of second-order 
curvature corrections and their application to the dark energy universe. Our starting 
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Figure 6. Variation of H, p, uj x and cu 2 with £ = +2, w = —1.5 and Q = 0. Initial 
conditions are chosen to be Hi = 0.2, fa = 2.0 and pi = 0.1. The Hubble rate continues 
to decrease while the energy density p increases. 
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Figure 7. Evolution of H and </> with £o = —2, w = —0.8 and Q — +5. Initial 
conditions are chosen to be Hi = 0.5, cfii = 2.0 and pi = 0.3. We find that both H and 
continue to decrease even in the presence of positive Q. 
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action (JTJ) is based upon the low-energy effective string theory with a barotropic 
perfect fluid coupled to a scalar field <fi. This is regarded as either the dilaton field 
characterising the strength of the string coupling, or a modulus field corresponding to 
the compactification radius of a dimensionally reduced theory. We adopted a general 
form of the second-order correction terms, given by equation (J2J), which includes the case 
of the Gauss-Bonnet curvature invariant. We have clarified the property of cosmological 
solutions in three classes of models: (i) fixed scalar field (section 0J), (ii) dilaton field 
(sections Eland EI), and (Hi) modulus field (section EJ. 

When the field <fi is fixed (model (i)), the existence of pure geometrical de Sitter 
solutions requires a spatial dimension d ^ 3, and that the parametrization of the second- 
order curvature correction is different from the GB one. Therefore these solutions are 
not realistic when applied to inflation or dark energy. Nevertheless we studied the 
classical stability of such solutions in the presence of a perfect fluid or a cosmological 
constant for completeness. 

If the dilaton field is not fixed (model (ii)), it is possible to construct pure dS 
solutions in the string frame provided that the evolution of the dilaton is linear in time 
(cf) = const). In the case of d = 3 spatial dimensions we obtained the dS solution given by 
Eq. ()60|) for the GB case, which was found to be stable. We also clarified the stability 
of other cases different from the GB parametrization. Moreover we showed that the 
existence of a phantom fluid generally makes the dS solution unstable. 

When a modulus field is present with fixed dilaton (model (Hi)), there exist 
cosmological solutions in which the field exhibits a logarithmic evolution, see equations 
(I70a|) - ()70e|) . We showed the existence of low-curvature solutions (|H3|) and Minkowski 
solutions 1)87)1. which can be joined each other if the coupling constant 5 given in 
equation (j24dj) is negative. In addition we obtained an exact solution (jl07|) for the 
modulus system, but this is found to be unstable. In the asymptotic future the solutions 
tend to approach the low-curvature one given by equation (|83jl rather than the others, 
irrespective of the sign of the modulus-to-curvature coupling 5. We also constructed 
high-curvature asymptotic solutions (jQB|) in the presence of a growing perfect fluid or 
a cosmological constant. Finally, we placed constraints on the viability of modulus- 
driven solutions using the current observational data. It is interesting that the GB 
parametrisation is excluded in any of the above mentioned regimes when a barotropic 
fluid is vanishing; see table [T] 

In sectionElwe solved the full equations of motion to simulate a dark energy universe 
in the future. The numerical results were then compared to the solutions obtained in 
previous sections. We used the modulus coupling with the GB curvature correction 
in 3 spatial dimensions. In the presence of the coupling Q between the field and 
the phantom fluid (w < —1), it is possible to consider situations in which the energy 
density of the fluid decays in the future. In fact we have numerically found that the 
big rip singularity can be avoided for the coupling Q which satisfies the condition 
Qu>2 — 3(1 + w)oj\ < —2 asymptotically (see figure EJ). This is actually achieved 
irrespective of the sign of Q and the asymptotic solutions are described by the low- 
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curvature one given by equation ()83j). When 5 is positive, we also found that even for 
Q = the Hubble parameter can decrease inspite of the growth of the energy density of 
the phantom fluid. Therefore the big rip can be avoided in this case as well, as illustrated 
in figure IH1 When the equation of state of the fluid is non-phantom (w > —1), we found 
that the presence of the coupling Q cannot give rise to the big rip singularity, contrary 
to what one might expect when considering the dark energy continuity equation (see 
figure Ej) . 

Another topic which is worth studying further is the behaviour of the asymptotic 
solutions at high redshift. For instance, when the deceleration parameter is parametrized 
as q(z) = qo + qiZ,§ supernovae observations give constraints for both go (redshift z ~ 0) 
and q\ (z ~ 0.5). Here we limited the discussion to the low-redshift evolution, since 
the experimental uncertainties on high-redshift data, together with the correlated data 
analyses, are more delicate to deal with. This topic is still under investigation. 

In the presence of a well-motivated underlying theory, one should specify the 
physical scale at which high-curvature terms with a non-GB parametrization arise and 
make sense, either at early times (high-energy inflation) or late times (dark energy 
around a big rip). This energy scale is related to the string coupling governing the 
perturbative higher-order corrections to the Einstein-Hilbert action. Although our main 
interest for the dark energy model of section 13 was only the Gauss-Bonnet case, we are 
still missing the precise theoretical setup in which a non-GB parametrization arise at 
high energies. This problem is still open and potentially rich of consequences. 

Also, in order to get close to real- world cosmology, we implemented a matter 
contribution into the Lagrangian and interpreted it as a dark energy source. Although 
the dilatonic/modulus corrections are motivated by string theory (in particular with 
the GB parametrization), the origin of the matter term still lacks in theoretical 
background. A completely coherent picture of dark energy should clarify its origin 
from first principles. Nonetheless our results do not rely on a particular choice for the 
kind of fluid, except for the assumption of the continuity equation (fTlj) and w = const. 
Field/string theory should also determine the coupling Q between the scalar field and 
matter. In our examples of big rip avoidance, the coupling Q is of order unity and has 
"natural" values anyway (in the sense specified by Dirac). 

We hope that this kind of phenomenology would be able to originate from a more 
precise formulation of the underlying physics. 
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The equations of motion ()33 j) - (j37jl . linearized under a linear perturbation of the vector 
X = (x, y, u, v, zY, read 5X = M5X, where the 5x5 matrix elements ra# are 

m u = mi 3 = m u = m 15 = 0, m 12 = l, (A.l) 
d — 1 ^ Ac\vx ^ 3(d — 3)cix 2 ze u 3a 4 v 4 y 2 ^ v 2 

2Xc2 C2 2c 2 d\c2X 2 2dc 2 x 2 2x 2 2d\c 2 x 2 ' 

y ze u 

m 22 = v - dx + -, m 23 = zm 25 = — , (A. 2) 

x d\c 2 x 

1 6a 4 t> 3 1c\x 2 v , . „ s 

m 24 = — + + — + y, (A.3) 

AC2 CXC2X C2 dAC2X 

m 31 = m 32 = m 33 = m 35 = 0, m 34 = 1, (A.4) 

d(d + l)x — dv + 2d(d + l)Acix 3 + MXcixy + 2dAa 4 t> 3 

m41 = T^6A^ ' (A ' 5) 

d(l + 2\cix 2 + \c 2 y) Qze u . . . 

m42 = 1 - 6Aa 4 , 2 ' m43 = " m45 = l-6Aa 4 t," (A ' 6) 

1 f , q 6Aa 4 u 

m 44 = < v — ax — 6Aa 4 f + 6aAa 4 f x ^ 2cm; 

1 — 6Aa 4 -ir [ 1 — 6Aa 4 t> 1 

— 2dy — v 2 — d(d + l)Acix 4 — d(d + \)x 2 — (4c 4 x 2 + c 2 y)d\y 

-4d\a 4 v 3 x + 3Xa 4 v 4 - Qze u ] J , (A.7) 

17151 = —d(l + w)z, m 52 = m 53 = 0, (A. 8) 

^54 = Q z , m 55 = —d(l + w)x + Qv, (A. 9) 

where we have dropped the subscript c for the background quantities. For the dS solution 
with Q — 0, two eigenvalues of the matrix M are given by 71 = and 72 = WI55. The 
other three can be extracted from 

7 3 - (m 2 2 + m 44 )7 2 + (m 22 m44 - m 24 m 42 - m 2 i)7 + (m 2 im 44 - m 24 m 4 i) = 0. (A. 10) 

In the four-dimensional case (d = 3), by inserting A = 1/4, a 4 = —1, c\ = 2, and the 
values found in equation (|dTI|). we find that the real part of the eigenvalues reads 

Re( 73 ) = 1.96, Re(7 4 ) = -2.40, Re(7 5 ) = -0.45, (A.ll) 

for C2 = 1, while 

Re( 73 ) = -0.45, Re( 74 ) = Re(7s) = -0.22, (A.12) 

for C2 = —1. We have checked that the sign of the eigenvalues does not change for 
different absolute values of c 2 . 

The GB case (c 2 = 0) must be treated separately. Together with the background 
equations (}3l)j) - (p!Tj) . we consider the (differentiated) Friedmann equation 

x = —v + (vz + z) — , (A.13) 
-£>i -t>i 
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where 

B x = d(d - l)x -dv + 2d{d - 3)AciX 3 - 6d\ Cl vx 2 , (A.14) 
B 2 = dx - v + QXa A v 3 + 2d\c l x 3 . (A.15) 

At the dS point the linear perturbations of the Bi coefficients are multiplied by vanishing 
factors. Defining the reduced vector X = (x, u, v, z)*, it can be shown that the 
eigenvalues of the 4x4 perturbation matrix M at the dS point are 

7i,3 = 0, 73 = m 55 = m 55 , (A.16) 

74 = mu + m 44 , (A. 17) 

where 

B2 _ B 2 m^\ 

™ n = ST 1 " 1 = Sx-lW (A ' 18) 
B 2 B 2 m M 

mi4 = BT AA = B x -B 2 m, 2 (A - 19) 

Note that m n m4 4 — m 41 m 14 = 0. If d — 3 and A = 1/4, the eigenvalue 74 is —0.45 for 

<24 = —1 and Ci = 2, and negative for other choices of parameters. Therefore the dS 

solution is always stable. Note that we may derive the same results by just multiplying 

the Jordan equation (|A.10|) times c 2 and then setting c 2 = 0. Since its coefficients are of 

the form (constant^ + (const ant ) 2 /c 2 , only the (constant) 2 terms survive and the above 

equation becomes linear. We have checked this shortcut explicitly. 

In the presence of a cosmological constant, z = A = const, there exist dS solutions 

only in the asymptotical limit m, t; -» as i - > 00 as explained in the text. The only 

difference with respect to the previous expressions for c 2 = is that they are evaluated 

at v — and rn.14 — > fhu + A/(Bi — B 2 m<i 2 ), 77144 — > rh^ = 77144 + Km^j [B\ — B 2 m^ 2 ). 

In this case the equation for the eigenvalue has a non-vanishing constant coefficient and 

the solutions are 

1 



73,4 



(mn + ^44) ± \l ("in + ^44) 2 + 4m A 

, ) 

Am 4 i 



(A.20) 



For negative Ci, one of the eigenvalues is always positive and the solution is unstable. 
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